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ABSTRACT 

Using reconstructed galaxy star formation histories, we calculate the instantaneous efficiency of galaxy star 
formation (i.e., the star formation rate divided by the baryon accretion rate) from z = 8 to the present day. This 
efficiency exhibits a clear peak near a characteristic halo mass of 10 1 L7 M , which coincides with longstanding 
theoretical predictions for the mass scale relevant to virial shock heating of accreted gas. Above the charac- 
teristic halo mass, the efficiency falls off as the mass to the minus four-thirds power; below the characteristic 
mass, the efficiency falls off at an average scaling of mass to the two-thirds power. By comparison, the shape 
and normalization of the efficiency change very little since z = 4. We show that a time-independent star forma- 
tion efficiency simply explains the shape of the cosmic star formation rate since z = 4 in terms of dark matter 
accretion rates. The rise in the cosmic star formation from early times until z = 2 is especially sensitive to 
galaxy formation efficiency. The mass dependence of the efficiency strongly limits where most star formation 
occurs, with the result that two-thirds of all star formation has occurred inside halos within a factor of three of 
the characteristic mass, a range that includes the mass of the Milky Way. 
Subject headings: dark matter — galaxies: abundances — galaxies: evolution 



1. INTRODUCTION 



Theorists have long predicted that galaxy formation is most 
efficient near a mass of ~ 1O 12 M based on analyses of su- 
pernova feedback, cooling times, and galaxy number counts 
(lSilk|1977||Rees & Osfriker|1977[|Dekel & Silk|1986l|White| 



& Rees|1978| |Blumenthal et al.|1984[ >. More recently, hydro- 
dynamical simulations have indicated that the host dark mat- 
ter halo mass strongly influences gas accretion onto galaxies 
UBirnboim & Dekel|[2003) |Keres et al.| |2005 1 De kel & BiriT| 
boim|2006 1. For low halo masses, these simulations predict 



that gas accretes in cold filaments ("cold mode accretion") 
directly to the galaxy disk, efficiently forming stars. Above 
a transition halo mass of 10 11 to 10 115 Mq (which is pre- 
dicted to be redshift-independent for z < 3), a shock develops 
at the virial radius which heats accreting gas ("hot mode ac- 
cretion") an d rapidly quenches star formation (Dekel & Birn- 
|boim|2006| ). 

To test these predictions, we use statistical reconstructions 
of the galaxy — halo connection for all observable galaxies 
( |Behroozi et al.||2012| l to compare the average star forma- 
tion rate in galaxies to the average baryon accretion rate as 
a function of halo mass and time. This approach allows us 
to directly test for a characteristic mass scale in the efficiency 
of star formation in halos. We summarize our methodology 
in present our main results in § |3"| and conclude in fflj 
Throughout this w ork, we assume a |Chabrier| ([2003 ) initial 
mass function, the |Bruzual & Charlot| ( |2003|l stellar popula- 
tion sy nthesis model, and the dust model in Blanton & Roweis 
( |2007[ ). We additionally assume a flat, ACDM cosmology 
with parameters £Im = 0.27, fl\ = 0.73, h = 0.7, n s = 0.95, 
and <j% = 0.82. 



2. METHODOLOGY 



In our main approach (fully detailed in Behroozi et al. 
|2012| l, we link galaxies observed at different redshifts to dark 
matter halos in a simulation so as to match galaxy stellar mass 
functions, average specific star formation rates of galaxies, 
and the cosmic star formation rate from z = to z = 8. Self- 
consistency is both important and difficult to achieve: ob- 
served star formation rates and known halo-halo merger rates 
constrain how galaxies placed into halos at one epoch can 
evolve into galaxies at a later epoch. We know, however, that a 
relationship for the average galaxy stellar mass at a given halo 
mass and redshift (SM(M,z)) must exist in the real Universe. 
We therefore adopt an extremely flexible (~ 20 parameter) 
fitting function for SM(M,z) so as to capture a wide range of 
possibilities for the true form and use an advanced Markov 
Chain Monte Carlo (MCMC) method to find the most likely 
functional forms in parameter space. As constraints, we cal- 
culate galaxy stellar mass functions (implied by SM(M,z) and 
the halo mass function) and galaxy star formation rates (im- 
plied by the time derivative of SM(M,z) and halo-halo merger 
rates) for each point in parameter space and compare them to 
observed results from z = to z = 8. We calculate uncertain- 
ties from a wide range of statistical and systematic effects (in- 
cluding uncertainties from stellar population synthesis mod- 
els, dust models, stellar population history models, sample 
variance, cosmology, photometric redshift errors, the faint- 
end slope of the stellar mass function, scatter between stella r 
mass and halo mass, etc.; see |Behroozi et aL |2010| |2012[ ), 
addressing concerns that our results may be biased by, e.g., 
limited observational constraints at high redshifts. 

We use observational co nstraints from 4 1 papers published 
over the last five years (see Behroozi et al. 2012| for a f ull list). 
These include res ults from SDSS and from PRIMUS ( |Mous-| 
|takas et aL]|2012| l, which uses low-resolution spectroscopy 
to accurately and consistently recover stellar mass functions 




FIG. 1 . — Top-left panel: Star formation rate as a function of halo mass and cosmic time in units of A/Q/yr. The grey shaded band excludes halos not expected 
to exist in the observable universe. Top-right panel: Conditional star formation rate as a function of halo mass and cosmic time, in units of the maximum star 
formation rate at a given time. Middle-left panel: baryonic mass accretion rate (MA) in halos as a function of halo mass and time, in units of M^lyr. Middle- 
right panel: the star formation rate to stellar mass ratio, in units of yr , as a function of halo mass and time. There is a roll-off towards higher halo masses; 
however, the normalization and characteristic mass are strongly redshift-dependent. Bottom panel: instantaneous star formation efficiency (star formation rate 
divided by baryonic mass accretion rate) as a function of halo mass and time. 
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FIG. 2. — The instantaneous star formation efficiency compared to the in- 
tegrated star formation efficiency at z = (i.e., stellar mass over the product 
of the baryon fraction with the halo mass). The shaded bands around each 
line show the one-sigma uncertainty contours. The integrated efficiency has 
a different peak and profile, as discussed in the text. 

from z = to z = 1 over a wide area of the sky. At high red- 
shifts, we include recent measurem ents of stellar masses and 
star formation rates to z = 8 (|Bouwens et al.|201 la|b||McLure| 



et al.||20PT Bradley et al. 2012 1. Notably, measurements of 



the cosmic star formation rate now agree with the evolution 
of the stellar mass density ( |Behroozi et al.|20T2}|Moster et al.| 
[20T2l|Bernardi et al.|2010p . 

For simulation data (halo mass functions, merger rates, and 
accretion histories), we extensively use the Bolshoi simulation 
( Klypin et al.|201 1 1. This simulation follows 2048 3 particles 
in a periodic, comoving volume 250 h~ l Mpc on a side using 



the art code ( |Kravtsov etal.|1997l|Kravtsov & Klypin|1999| ); 
it has a mass resolution of 1.9 x 10 Mp, and a force resolution 



of 1 h~ l kpc. The adopted flat, ACDM cosmology (see 
is consistent with the latest WMAP7+BAO+H results ~ 
|matsu et al. |20 1 0) > . Simula tion analysis was perfo rmed using 
the ROCKS TAR halo finder (jBehroozi et al.|201 la[ > and merger 
tree code in Behroozi et al.| ( |201 lt>J >. 

3. RESULTS 

3.1. Strong Mass Dependence for the Star Formation 
Efficiency 

We show the main output of the method in |Behroozi et al.| 
(2012), the average star formation rate (SFR) in dark matter 
halos as a function of virial mass (Bryan & Norma njl998| l and 
time, in the top-left panel of Fig.[T] The SFR depends strongly 
on time, yet there is also a distinct embedded mass scale, as 
may be seen by normalizing to the maximum SFR as a func- 
tion of time (Fig. |T| top-right panel). To understand the impli- 
cations for gas physics in halos, it is necessary to consider the 
baryon accretion rate as well. Thi s rate is calcula ted here as 
the halo mass accretion rate (Behr oozi et al.|2012[) ti mes the 
cosmic baryon fraction (see |van de Voort et al.|201 1} , and is 
shown in Fig.[T] middle-left panel. 

The baryon accretion rate increases with halo mass and 
lookback time. This trend combines with trends in the star 
formation rate to reveal a clear picture of star formation ef- 
ficiency in halos (Fig. [T] bottom panel). This efficiency, de- 
fined as the star formation rate divided by the baryon accre- 



tion rate, shows a prominent maximum near a characteristic 
mass of 10 1L7 M Q (see also Fig. [2). Indeed, the star forma- 
tion efficiency over 90% of the history of the universe (z < 4) 
is strongly dependent on halo mass; by comparison, it has a 
weak dependence on time. 

The peak in the star formation efficiency (SFE) at 10 1 1 7 M 
represents observationally-constrained evidence for a charac- 
teristic mass for galaxy formation. This characteristic mass 
matches longstanding theoretical predictions in both its value 
and in its lack of evolution since z = 3 — 4. The steep effi- 
ciency cutoff above the characteristic mass (SFE cx M h 4 ^, 
where M/, is halo mass) suggests that a strong physical mech- 
anism prevents incoming gas from reaching galaxies in mas- 
sive halos. Besides the effect of hot-mode accretion, this slope 
coincides with the mass and luminosity scaling for supermas- 
sive black holes (Lbh <x Mbh', scalings f or Mbh vary from 



M BH oc a 4 oc M^ 3 to M BH oc a 5 <x M 5 J 3 ; |Ferrarese & Mer- 



|ritt|2000[ pVIcCo nnell et a l.|201 l| l, which may prevent residual 
cooling flows in massive clusters from forming stars. Below 
the characteristic mass, the efficiency is not a perfect power 
law; between Mf, ~ 10 10 and ~ 1O 1L5 M , the average slope 

is SFE oc M 1 ^ . This may seem to be consistent with semi- 
analytic galaxy formation models that use supernova feedback 

(most commonly scaling as V} lrr o c M^ 3 ; |Luetal.|2012[|Hat-| 
|ton et al.]|2003| |Somerville et al. 2008) to expel most or. trie 
gas in low-mass halos. However, these models often assume 
that the expelled gas rea ccretes onto the halo after a dynami- 
cal time ( Luetal.|2011| l; this extra incoming gas would result 
in a steeper mass-dependence for the SFE at low halo masses. 

3.2. Weak Time Dependence of the Star Formation Efficiency 

The weak time dependence of the star formation efficiency 
is unexpected given the different environments of 10 11 7 M Q 
halos at z = 4 and z = 0. At z = 4, the background matter den- 
sity was ~ 125 times higher, mass accretion rates were ~ 40 
times higher, galaxy-galaxy merger rates were ~ 20 times 
higher, and the UV background from star formati on was ~ 
500 times m ore intense than at the present day (Behroozi 
et al. l 2012J I. None of these differences significantly influ- 
enced star formation efficiency (unless they conspired to can- 
cel each other out), strongly constraining possible physical 
mechanisms for star formation in galaxies and halos. 

While the time dependence of the SFE is weak, it is not 
absent. As seen in Fig.[T] the characteristic halo mass evolves 
from a peak of 10 12 M Q at z = 3 to 10 1 1 5 M Q at z = 0. The peak 
star formation efficiency also evolves around its average value 
of 0.35, reaching a maximum of 0.55 at z = 0.8 and a minimum 
of 0.22 at z = 0. However, observational constraints on star 
forma tion rates and stellar masses are uncertain at the 0.3 dex 
level ( |Behroozi et al.|2012| |2010| i especially for z > 1 ; these 
are larger than the observed deviations (±0.2 dex) in the peak 
star formation efficiency. The variations in the characteristic 
mass are likely more significa nt; while observational biases 
can be stellar mass-dependent (Behroozi et al.|20T2| in a way 
that chang es the location of the peak halo mass, this effect 
(<0.1 dex |Leauthaud et aT7||2012| l cannot account for the 0.5 
dex change from z = 3 to z = 0. Nonetheless, these concerns 
do not alter the fact that the trends with mass (four decades 
of variation) in the star formation efficiency are stronger than 
the trends with time. 

One way to eliminate the residual time dependence in the 
characteristic mass is to use a different mass definition. For 
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FIG. 3. — Left panel: the stellar mass to halo mass ratio at multiple redshifts as derived fr om observa tions (Behroozi et al. 2012) compared to a model which 
has a time-independent star formation efficiency (SFE). Error bars show 1 — a uncertainties (Beh roozi et al.|2U12^ . A time-ind ependent SFE predicts a roughly 
time-independent stellar mass to halo mass relationship. Right: the cosmic star formation rate tor a compilation or observations ( Behroozi et al. 2012 1 compared 
to the best-fit model from a star formation history reconstruction technique (Behroozi et al. 20121 as well as the time-independent StH model, the latter model 
works surprisingly well up to redshifts of z ~ 4. However, a model which has a constant efficiency (with mass and time) also reproduces the decline in star 
formation well since z ~ 2. 
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FIG. 4. — Left panel: Star formation rate as a function of halo mass and cosmic time, weighted by the number density of dark matter halos at that time. Contours 
show where 50 and 90% of all stars were formed; dashed line shows the median halo mass for star formation as a function of time. Right panel: Star formation 
rate as a function of galaxy stellar mass and time, weighted by the number density of galaxies at that time. Contours and dashed line are as in top-left panel; 
dotted line shows current minimum stellar masses reached by observations. 



example, using Mxoob (i-e-, 200 times the background density) 
would cancel some of the evolution from z = 1 to z = 0. How- 
ever, this would also raise the mass accretion rate at z = 0, 
which would increase evolution in the star formation effi- 
ciency's normalization. Using the maximum circular velocity 
(Vcirc) or the velocity dispersion (a) instead would also lead 
to more evolution in the SFE (at fixed V c j re or a): due to the 
smaller physical dimensions of the universe at early times, 
both these velocities increase with redshift at fixed virial halo 
mass. 

The nearly-constant characteristic mass scale is robust to 
our main assumption that the baryon accretion rate is propor- 
tional to the halo mass accretion rate, because this mass scale 
is already present in the conditional SFR (Fig. [TJ. A baryon 



accretion rate which scales nonlinearly with the dark matter 
accretion rate would change the width of the most efficient 
halo mass range, but it would not change the location. How- 
ever, as discussed previously, the baryon accretion rate for 
small halos (M h < 1O 12 M ) can differ from the dark matter 
accretion rate through recooling of ejected gas; the changing 
virial density threshold can also introduce non-physical evolu- 
tion in the halo mass which affects the accretion rate ( |Diemer| 
et al. 2012| l. Properly accounting for these effects may change 



the low-mass slope of the star formation efficiency; we will 
investigate this in future work. 

Note that the level of consistency seen in the star forma- 
tion efficiency is not possible to achieve using other common 
specific ratios, e.g., the specific star formation rate (SFR to 
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stellar mass ratio; Fig. [T] middle -right panel) or the SFR to 
halo mass ratio. The stellar mass — halo mass ratio (i.e., the 
integrated form ation efficiency) does show somewhat simi- 
lar features (Leauthaud et al.|2012} Moster et al.|2012||Wang| 
eTaT1|20T2| IBehroozi et al |)2012| ]2T)10| |Yang et aL||2012[ 



Conroy & Wechsler 2009 )Thowever, the integrated efficiency 
is several steps removed from the actual physics of star for- 
mation. Galaxy stellar mass is influenced by stellar death, 
galaxy-galaxy mergers, a nd ejection of merging stellar mass 
into the intracluster light (Behroozi et al.|2012| Moster et al.| 
|2012| |Conroy & Wec hsler 2009J), complicating the interpre- 
tation of the integrated efficiency. Moreover, the shape of the 
integrated efficiency is influenced by star formation along the 
entire halo mass accretion history. Thus the profile of the in- 
tegrated efficiency differs notably from the instantaneous star 
formation efficiency which we study here, shown in Fig. [2j 
especially in terms of the peak efficiency and the high-mass 
slope. 

3.3. A Time-Independent Model 

Going further, it is interesting to approximate the star for- 
mation efficiency for individual halos as completely time- 
independent. In this case, the stellar mass at a given halo mass 

r'dM h „dSM, f M "dSM iA/r 

-fb-j7T dt = fb I -i^rd M h (1) 



SM = 



dt dM h 



dM h 



(where SM is the stellar mass, M/,(t) is the halo mass accretion 
history, and fb^f is the star formation efficiency). The total 
stellar mass formed then becomes a function of only the halo 
mass and not of time. 

The specific choice of redshift for the instantaneous star for- 
mation rate does not matter greatly, as shown in Fig. [T] We 
nonetheless marginalize the instantaneous star formation rate 
over time; the resulting functional form is shown in Fig. [2] 
Using this as the time-independent efficiency, we calculate 
the total stars formed as a function of halo mass using Eq. [T] 
and reduce the resulting value by 50%, corresponding to the 
stellar population rem aining for a 6 Gyr-old starburst ( |Con-| 
roy & Wechsler 2009). (For comparison, a 1 Gyr-old starburst 
would have 60% of its original stars remaining). This allows 
us to calculate the stellar mass to halo mass ratio, as shown in 
the left panel of Fig. [3] Similarly, we may use halo mass ac- 
cretion rates and number densities along with the same SFE to 
calculate the cosmic star formation rate (Fig. [3] right panel). 

The real universe is more complicated, of course; the stel- 



curately reproduce galaxy number counts ([Leauthaud et al. 


2012 


Moster et al. 2012 


Wang et al.||2012| [Behroozi et al. 


2012 


2010; Yang et al. 


2012 Conroy & Wechsler| |2009 ; 


Moster et al. 2010). However, integrating a time-independent 



SFE with respect to halo mass reproduces the z = stellar 
mass to halo mass relation to within observational systemat- 
ics over nearly five decades in halo mass (10 10 to 10 15 M Q ). 
Similarly, integrating the SFE times the mass accretion rate 
and number density of halos gives a precise match to the ob- 
served cosmic star formation rate from z ~ 4 to the present. 

Furthermore, the prediction in time-independent SFE mod- 
els of fixed stellar mass formed at a given halo mass is not 
far off from observational constraints at z = (0.2 dex scatter 
in stellar mass at fixed halo mass; Red dick et al.|20i2| ). The 
evolution in the median stellar mass to halo mass relation with 
time, corresponding to an evolution in the star formation ef- 
ficiency, may then set a lower bound on the scatter in stellar 



mass at fixed halo mass at the present day. Conversely, the 
scatter in stellar mass at fixed halo mass today sets an upper 
bound on the possible evolution of the median stellar mass to 
halo mass ratios at earlier times. 

As regards the cosmic star formation rate, it is tempting to 
view the time-independent efficiency model as more of a suc- 
cess (in terms of matching galaxy formation physics) than it 
actually is. In fact, a model in which the star formation effi- 
ciency is 7% regardless of halo mass or time can also match 
the observed decline in cosmic star formation rates (Fig. [3] 
right panel), but would not match the stellar mass to halo mass 
ratio or galaxy number counts. For that reason, the decline in 
the cosmic star formation rate since z = 2 is more related to de- 
clining dark matter accretion rates than changes in how galax- 
ies form stars. This may explain past successes in reproducing 
the cosmic SFR with a variety of in compatible physica l mod- 
els (e .g., Krumholz & Dekel|2012| |Dave et al.|2012|[Bouche] 
|et al.p0~l0| |Hernquist & Springel 2003 ). That said, the rise 
in the cosmic star formation rate from early times to z = 2 is 
much steeper than a mass-independent efficiency model pre- 
dicts. Matching this rise is much more closely tied to galaxy 
formation physics, as it requires an increase in the average star 
formation efficiency with time. In the mass-dependent model, 
this is provided by an increasing number of halos reaching the 
characteristic mass. 

3.4. Consequences for When and Where Stars Were Formed 

The star formation efficiency leaves a distinct imprint on the 
star formation history of the universe: as halos pass through 
the characteristic mass (1O 1L7 M ), they form most of their 
stars. Equivalently, most stars were formed in halos between 
10 n 5 M Q and 10 122 M© (Fig.H left panel). Furthermore, be- 
cause of the tight correlation between stellar mass and halo 
mass, most stars formed in galaxies with stellar masses be- 
tween 10" and 1O 1O 8 M (Fig. right panel). This same 
narrow range of halo and stellar masses ( which includes the 
stellar and halo mass of the Milky W ay; Busha et al.||201 1 



Flynn et al.|2006| |Klypin et al.|2002} is responsible for most 
star formation since at least z = 4, due to the constancy of the 
star formation efficiency with time. Given current observa- 
tional limits (Fig. [4] right panel), surveys have probed a stel- 
lar mass and redshift range corresponding to 90% of the star 
formation in the Universe. 

4. CONCLUSIONS 

As we have shown, the efficiency of galaxy formation falls 
off strongly on either side of a characteristic halo mass and, 
in the absence of a conspiracy of effects, is weakly correlated 
with time and environment. This would suggest a model for 
galaxy formation in which self-regulation after z ~ 4 is nearly 
perfectly efficient and is controlled by effects which corre- 
late largel y with the loca l gravitational potential: supernova 
feedback (Dekel & Silk 1986 ) and possibly metallicity effects 
( jKrumhoIz & Dekel||2012| )Timit galaxy growth in low-mass 
halos, and hot mode accretion as well as black hole feedback 
(Silk & Rees 1998) limit growth in high-mass halos. Under- 
standing how these physical models for feedback are incom- 
plete, especially in terms of the steep efficiency falloff on ei- 
ther side of the characteristic mass and its mild evolution with 
time, will remain a challenge for future research. 
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